The heat flux from a relativistic kinetic equation 
with a simplified collision kernel 
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Abstract. We show how using a special relativistic kinetic equation with a BGK- like collision 
operator the ensuing expression for the heat flux can be casted in the form required by Classical 
Irreversible Thermodynamics. Indeed, it is linearly related to the temperature and number density 
gradients and not to the acceleration as the so-called "first order in the gradients theories" contend. 
Here we calculate explicitly the ensuing transport coefficients and compare them with the results 
obtained by other authors. 
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The purpose of this work is to calculate the explicit form for the heat flux in an ideal 
relativistic gas to find out what is the form of the constitutive equation as required by 
Classical Irreversible Thermodynamics (CIT), in the temperature and particle number 
density representation. This result is important by itself because it will help in clarifying 
whether or not such constitutive equation should contain a term proportional to the 
acceleration as it has been claimed by many authors. To accomplish this task, we start 
by considering the simplified relativistic Boltzmann equation in the absence of external 
forces: 

The term on the right hand side of Eq 
by means of the BGK approximation yj] 
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where w £ is the molecular three-velocity and y = (l — w l Wi/c 2 ) 1/7 is the usual rel- 
ativistic factor. All greek indices run from 1 to 4 and the latin ones run up to 3. The 
number density reads 

n = J fydv* (3) 



J(f,f) 



(1) 



. (OQ) is the collision operator which is modeled 
. In Eq. (OQ), the molecular four velocity, v a is 
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with dv* = y 5 10]. Now, standard kinetic theory leads to the energy balance equation 
Oil: 
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where the thermodynamic average is defined through 
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In Eq. © one can easily identify the heat flux as: 

jf Q] = n (mc 2 v l } = mc 2 J v l fydv* (6) 

In a Chapman-Enskog expansion, the distribution function is written as 

/ = /(°)(l + 0) (7) 

where is the first order correction in the Knudsen parameter, a weighted measure of 
the gradients in the system. For particles of rest mass m, relativistic parameter z = £fz 
and in the non-degenerate case, the equililibrium function reads OH: 

/ (0) = / (8) 

in which k is Boltzmann's constant, the hydrodynamic four-velocity and Jfa is the 
modified Bessel function of the second kind. We remind the reader that the first order 
correction contains already the dissipative effects in the system. 

If T denotes the collisional time appearing in the BGK model, using the functional 
hypothesis [4J], the function in Eq. © is given in terms of the thermodynamic forces 
as 

t = -rv \—n ta + - sr T, a + -j- jr1 &) (9) 



where the derivatives of the Jiittner distribution function are 
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FIGURE 1. The functions f T (z) (left) and /„ (z) (right). 



We now introduce Euler's equations for the time derivatives in ©. A hydrodynamic 
acceleration term appears in the last term of Eq. © for a = 4. In the standard approach 
this term is expressed by the pressure gradient using Euler's equation. This would lead to 
an awkward representation involving n, T, p and u a . However, in order to comply with 
the requirements of CIT, using the fact that the ordinary non-equilibrium state variables 
of the fluid are n, u a and T, we substitute the pressure gradient through the equation of 
state p = nkT, which readily leads to the following form for the heat flux 
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J [Q] = ~ Ltt ^ +L " t ~ (13) 

Equation (fT3l) is of the canonical form required by CIT, the thermodynamic forces are 
the temperature and the density gradients, respectively. The transport coefficients are 
given by 

nk 2 T 2 , . nk 2 T 2 

Lt = Tfr (z) , Ln = Xfn z (14) 
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where the functions fj and are given by 

«>='^*f§)lMi 



! + ^7TT+ \\\ (15) 

and are shown in Figure 1 . 

The coefficient given in Eq. (fT6l) has not been analyzed before, although it has been 
proposed in Ref. |5|]. Its presence in the constitutive equation for the heat flux given in 
Eq. (fT3l) . which has the structure required by CIT, guarantees the validity of Onsager's 



regression of fluctuations hypothesis. Other authors obtain an expression for the heat 
flux in terms of the acceleration or pressure gradient which is, not only against the tenets 
of CIT, but also gives rise to the so-called generic instabilities for the relativistic fluid 
[0, 0]. In contrast, our approach provides a picture in which fluctuations do not grow 
exponentially within the first order in the gradients description of the relativistic simple 
fluid. 

Two important facts shown in Figure 1 are worth pointing out as a final remark. 
Firstly, we have obtained an analytical expression for the thermal conductivity for the 
relativistic gas, using in the BGK approximation. This result is new and differs from 
previous calculations since this is the first time the representation given in Eq. (fT3l) is 
explicitly used. Notice also how this coefficient yields the correct non-relativistic limit 

Lt ~ § r T. Secondly, the figure on the right shows a novel result, namely a coefficient 
relating a heat flux to a number density gradient present only in the relativistic case. This 
coefficient vanishes for z — > 0. 
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